Mon. Not. R. Astron. Soc. 000, [HS] (2011) Printed 1 March 2013 (MN I^T^ style file v2.2) 



TeV blazar variability: the firehose instability? 
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ABSTRACT 

Recently observed minute timescale variability of blazar emission at TeV energies has 
imposed severe constraints on jet models and TeV emission mechanisms. We focus 
on a robust jet instability to explain this variability. As a consequence of the bulk 
outflow of the jet plasma, the pressure is likely to be anisotropic, with the parallel 
pressure P|| in the forward jet direction exceeding the perpendicular pressure Pj_. 
Under these circumstances, the jet is susceptible to the firehose instability, which can 
cause disruptions in the large scale jet structure and result in variability of the observed 
radiation. For a realistics range of parameters, we find that the growth timescale of 
the firehose instability is « a few minutes, in good agreement with the observed TeV 
variability timescales for Mrk 501 (Albert et al. 2007) and PKS 2155-304 (Aharonian 
et al. 2007). 
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1 INTRODUCTION 

Several blazars have recently displayed minute-timescale 
variability at TeV energies (Albert et al. 2007; Aharonian 
et al. 2007). This has resulted in a great deal of interest in 
the emission mechanisms responsible for these TeV flares. 
Suggestions range from a coherent instability in a compact 
emission region (e.g., Begelman, Fabian & Rees 2008), mis- 
aligned minijets inside the main jet (e.g., Giannios, Uzden- 
sky & Begelman 2010), jet deceleration (Georganoupolos & 
Kazanas 2003; Levinson 2007) to wiggles in an anisotropic 
electron beam directed along the jet (Ghisellini et al. 2009). 
Considerable attention has also been paid to correlated vari- 
ability at other wavelengths (such as optical wavelengths) 
exhibited by TeV blazars (e.g., Gopal-Krishna et al. 2011 
and references therein). In general, there have been several 
mechanisms proposed for producing the observed variabil- 
ity in the jet emission, ranging from plasma mechanisms 
(Krishan & Wiita 1994) to beamed radiation (e.g., Crusius- 
Watzel & Lesch 1998). TeV variablity is somewhat unique, 
in that explaining even steady-state TeV emission from 
blazars poses significant challenges. Electron emission (via 
the synchroton self-compton process or the external inverse 
Compton process) is frequently invoked. While electrons are 
certainly the readiest radiators, both these mechanisms re- 
quire very highly relativistic electrons (random Lorentz fac- 
tors 7 ~ 10'* — 10^ or more). This in turn poses severe elec- 
tron reacceleration problems. 
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2 CURRENT MODELS FOR TEV FLARES: 
THE ROLE OF BULK OUTFLOW 

We now pay attention to the role of bulk plasma outflow 
in a couple of popular models for TeV flares. In the first 
model we discuss, the bulk plasma outflow is necessary to 
alleviate the problem of copius pair production (and con- 
sequent degradation of TeV photons) in the TeV emission 
region. In the second one, the idea of a highly anisotropic, 
directed beam of TeV emitting electrons is central to ex- 
plaining the observed variability. We point out later that 
the bulk Lorentz factor (F) of the plasma outflow is related 
to pressure anisotropy in the jet. The magnitude of F (and 
therefore that of the pressure anisotropy) is limited by the 
excitation of hydromagnetic waves. We will then examine 
the role of the pressure anisotropy in exciting the well stud- 
ied firehose instability in the jet. Jet disruption due to this 
instability is found to be a viable explanation for the ob- 
served variability at TeV energies. 



2.1 TeV variability due to a compact emission 
region: the pair production problem 

The first scenario we examine is due to Begelman, Fabian 
& Rees (2008). In this scenario, the TeV photons are pro- 
duced via inverse Compton upscattering of an existing soft 
photon population by energetic electrons. Using a peak fre- 
quency of lO^'' Hz for the soft photons and a bulk Lorentz 
factor of F ~ 50, they estimate that the random Lorentz 
factor of the energetic electrons 7 « 10* in order to produce 
TeV emission. They relate the short variability timescales 
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(of the order of a few minutes) to the size of the emission 
region via the hght crossing time argument. This argument 
implies that the emission originates from a region whose size 
is a small fraction of the black hole's Schwarzschild radius 
(Begelman, Fabian & Rees 2008). For relativistic, Poynting 
flux dominated jets, it is possible that structures as small 
as the gravitational radius of the black hole are imprinted 
on the jet as it is launched, and modulate the emission far 
away from the central region (Kirk & Mochol 2010). An- 
other possibility is that the TeV flares could originate in 
highly localized fluctuations in the outflowing jet plasma 
due to instabilities. 

Importantly, Begelman, Fabian & Rees (2008) show 
that the outflowing plasma must have a substantial bulk 
Lorentz factor F ~ 50 in order for the TeV photons to 
escape from the compact emission region without produc- 
ing pairs. If this was not so, copius pair production due to 
the interaction between the TeV photons and locally pro- 
duced synchrotron photons would substantially degrade the 
TeV photons produced. Furthermore, Begelman, Fabian & 
Rees (2008) argue that the TeV emitting region needs to 
be located far away from the central black hole (at least 
100 Schwarzschild radii) in order to avoid rapid degradation 
of the TeV photons due to pair production on the soft pho- 
ton background of the accretion disk (e.g., Becker & Kafatos 
1995). Assuming a reasonable jet opening angle, this implies 
that the size of the emission region is around two orders of 
magnitude smaller than the transverse dimension of the jet, 
at the distance where the variable TeV emission originates. 
This scenario thus envisages a very compact emission re- 
gion comprising energetic electrons flowing out with a bulk 
Lorentz factor « 50. 



2.2 TeV variability due to beamed electron 
distribution? 

Ghisellini et al. (2009) outline an interesting scenario where 
the TeV emission is due to a highly anisotropic electron 
beam. The electrons are almost co-aligned in a narrow beam. 
Another way of stating this could be to say that the electrons 
are hot along the beaming direction and cold in the perpen- 
dicular direction; the jet pressure (as observed by a distant 
'lab" observer) is highly anisotropic. It would be possible to 
define a meaningful bulk Lorentz factor along the beaming 
direction, but the jet pressure would be anisotropic even in 
the bulk comoving frame. The electrons could be emitting 
TeV radiation via the inverse Compton scattering of soft 
photons. As an extreme case, they calculate the electron 
Lorentz factor in a situation where there is no bulk motion 
and the soft photons are isotropically distributed (although, 
as mentioned earlier, the highly anisotropic electron beam 
clearly implies ordered bulk motion). They find that the 
electron random Lorentz factor 7 needs to be « 10^ in order 
to account for TeV photons produced via inverse Compton 
scattering off soft photons of energies ~ 1 eV. Ghisellini et 
al. (2009) envisage large-scale wiggles in the beam, so that 
the radiation would be observed only for the fraction of time 
that the beam points toward the observer, leading to vari- 
ability. The electron beam is assumed to retain its coherence, 
and the variability is only because of the fact that the large 
scale magnetic field points toward the observer for a limited 



time. Such a scenario avoids explicit reference to a compact 
emission region. 

A central feature of Ghisellini et al. (2009) 's model is 
the highly beamed, anisotropic electron distribution, which 
leads to the parallel pressure greatly exceeding the per- 
pendicular one. We show below that the electron pressure 
anisotropy is unlikely to assume extreme values. We show, 
however, that even modest amounts of pressure anisotropy 
are enough to destabilize the large-scale structure of the jet 
via the firehose instability, and this could result in observed 
variability. 



3 OUR SCENARIO 

Our scenario is similar to Begelman et al. (2008) 's one in 
that the radiating electrons possess a random Lorentz factor 
7 ranging from lO'' to lO'' as well as a bulk Lorentz factor 
r. The issue of pressure anisotropy highlighted by Ghisellini 
et al. (2009) is important in our scenario too. 

In this work, we envisage a jet that comprises an 
electron-proton plasma which streams along the forward 
jet direction. All the quantities referred to from here on- 
wards are defined in the distant observer's frame of refer- 
ence. Electrons would be bound to the protons, and stream 
along with them, but would also possess random Lorentz 
factors 7 ~ 10*-10^ . Since the electrons have such high ran- 
dom Lorentz factors, and the protons are most likely much 
colder, the electrons are the dominant contributors to the 
pressure. The velocity with which the plasma streams along 
the jet is typically restricted to the Alfven speed va via a 
self-limiting process. This limit on the streaming speed is 
imposed by hydromagnetic waves emitted by the electrons 
(Wentzel 1969; Melrose 1970), which also scatter the elec- 
trons and limit the degree of anisotropy. Following Wentzel 
(1968; 1969) we envisage the equilibrium particle distribu- 
tion function in momentum p of the jet plasma to be of the 
form 



/(p) = /.(p) 1 + 3!m 



,VA 



(1) 



where /i(p) is the isotropic part of the distribution function 
and satisfies 



I Mp)d'p = i. 



(2) 



The quantity jj, is the cosine of the particle pitch angle and 
represents the fraction of the particle momentum along the 
direction of the bulk motion of the jet. The quantity c is the 
speed of light. The second term in the brackets expresses 
the anisotropy of particle distribution. It may be noted that 
the anisotropy in this formulation is typically only a pertur- 
bation; in other words, we require that va/c <C 1. Wentzel's 
(1968) formulation includes a higher order term oc /i'^ on the 
right hand side of Eq[T] but we don't consider such a term in 
the interest of simplicity. While we do not appeal to a spe- 
cific mechanism to generate the electron anisotropy, such 
as the pitch angle dependence of synchrotron losses (e.g., 
Melrose 1970; Wentzel 1969) or the intrinsic jet launching 
mechanism (e.g., Aharonian, Timokhin & Plyasheshnikov 
2002) , we emphasize that this distribution (Eq [ij is only 
mildly anisotropic, and the degree of anisotropy in is re- 
stricted to « va/c ^ 1. Furthermore, the bulk streaming 
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speed of the jet is restricted to be approximately equal to 
the Alfven speed va- The ratio of the parallel pressure P| 
to the perpendicular pressure Pi_ in the jet arising from the 
distribution given by Eq[T]is (Wentzel 1968) 



PU_ ^ 1 + va/c 
P±^ I- va/c 



(3) 



Since the bulk streaming speed of the jet is ~ ua , we 
can write the following expression for the jet bulk Lorentz 
factor F: 



r = (i 



2 / 2n 

Va/c ) 



^1/2 



(4) 



keeping in mind the fact that F cannot greatly exceed unity, 
since va/c -ti 1. 

Before we move on to analyzing the firehose instability 
arising out of the anisotropic pressure distribution in the 
jet, we note that electron reacceleration within the jet is 
essential, in view of the copious radiative losses they expe- 
rience. The electrons could be energized by resonating with 
the very hydromagnetic waves they shed, much in the same 
manner as cosmic rays are thought to be accelerated (e.g., 
Bell & Lucek 2000), or by wave-particle interactions with a 
separate, pre-existing wave population (e.g., Eilek 1979). It 
is also possible that electrons are reaccelerated by the turbu- 
lence initiated by a previous episode of the firehose instabil- 
ity; such a scenario is often invoked for electron acceleration 
in solar flares (e.g., Paesold & Benz 1999). The protons, on 
the other hand, are likely to remain substantially colder. 



3.1 The Firehose instability due to Py > P± 

The well known firehose instability (also often referred to 
as the gardenhose instability) is typically operative when 
the parallel pressure P|| in the forward jet direction exceeds 
the perpendicular pressure P± transverse to the jet (e.g., 
Krall & Trivelpiece 1973). This nonresonant, fluid instabil- 
ity will result in large-scale displacements in the beam, and 
possible disruption of the large-scale magnetic field in the 
jet (e.g.. Baker et al. 1988), which could result in variabil- 
ity in the observed TeV emission. Resonant instabilties can 
also be operative, but they typically only result in a redis- 
tribution of electron pitch angles, while nonresonant, fluid 
instabilities disrupt the large-scale magnetic field. The argu- 
ments of the previous section suggest that the jet pressure 
is unlikely to be highly anisotropic. We show hereafter that 
the growth timescales for the firehose instability in a typical 
TeV emitting jet are in agreement with observed variability 
timescales, even for moderate pressure anisotropics. 



3.2 Growth timescale of the firehose instability 

While there are several treatments of the firehose instability, 
we concentrate on one that offers a convenient analytical ap- 
proximation to the maximum growth rate of the relativistic 
firehose instability from a solution of the dispersion equa- 
tion. The maximum growth rate lj^i of this instability can 
be approximated to within ~ 20% by the following expres- 
sion: (Noerdlinger & Yui 1969) 



= 0.55 WLR 



(1-g) 

8-f Ai/2 



(5) 



The quantity ulr is the mean relativistic gyrofrequency. 
given by 



wlr = 



eB 



(6) 



where e is the electron charge, B is the large-scale mag- 
netic field in Gauss, 7 is the random Lorentz factor of the 
TeV emitting electrons, rUe is the electron mass and c is the 
speed of light. The growth timescale given by Eq ([5} is an 
approximation to the growth rate derived from solving the 
dispersion relation for the instability. 

The quantity H is the magnetic pressure nondimension- 
alized by the particle anisotropy, and is defined as 



H ^ 



B^ 



47rP||(l-Px/P|| 
The quantity A is defined as 

A = l + -^, 



(7) 



(8) 



and can be related to the bulk Lorentz factor F via Eq Q. It 
may be noted that the gyrofrequency is used in Eq ([5]) only 
as a convenient parametrization for the growth timescale, 
and does not suggest any kind of resonance; the relativistic 
firehose instability is a nonresonant one. 
Identifying the Alfven speed as 



2 
Va 



B^ 



(9) 



where nip is the proton mass and TV is the particle num- 
ber density, we can express H (Eq[7| as 



H 



vl 



7C2(1-Px/P||) 



(10) 



where m^ denotes the electron mass. We have used the pro- 
ton mass in computing the Alfven velocity, since they are 
the heavier species, and are the primary contributors to the 
matter density. We have also used P| | — ■y N rrie c? in deriv- 
ing Eq (|10|l . since we expect the energetic electrons (with 
7 « 10'') to be the primary contributors to the jet pressure. 
The maximum growth timescale of the relativistic fire- 
hose instability is 



igr ^ 27r/a;g 



(11) 



As mentioned in § 3, the growth timescale is evaluated in 
the frame of the distant observer. Using Eqs (|10p . (|8l), (|6|, 
(O, dH) and @ in Eq HIJ, we get 



27r 



jmeC 



^' "" 0.55 eP 

1 TTl) 

7 "^e 



;+ ( i + (i-F"") 

l-F(l-F-2)i/^ 



1/2' 



2(l-F-2)i/2 



(12) 



4 Suhramanian, Shukla and Becker 




0.01 



0.02 0.03 

B (Gauss) 



0.04 



0.05 



Figure 1. The predicted variability timescale tgr in minutes as a 
function of the magnetic field B in Gauss. The thin solid line is 
for 7 = 5 X lO"*, r = 1.01 and the thick solid line is for 7 = 5 X 10*, 
r = 5. The thin dashed line is for 7 = lO", T = 1.01 and the thick 
dashed line is for 7 = lO'', F = 5. 



ratio of the Alfven speed to the speed of hght va/c <C 1, 
the bulk Lorentz factor F needs to be restricted to values 
that are not appreciably greater than unity (EqU). These 
values for F imply relatively modest constraints on the jet 
production mechanism. Although rather high values for F 
are often discussed in the literature, it is worth mentioning 
that multi-epoch radio mapping of parsec-scale jets in TeV 
emitting blazars suggest that F « 1-3 (Finer, Pant & Ed- 
wards 2008; also see Giroletti et al. 2004). Similar values 
for F are also implied by blazar unification schemes (e.g., 
Urry & Padovani 1991). Theoretical considerations regard- 
ing launching electron-proton jets imply that asymptotic 
bulk Lorentz factors are likely to be limited to values well 
below 10 (Subramanian, Becker & Kazanas 1999). 

The broad conclusion from Figure 1 is that a wide range 
of realistic parameters (5 x 10* < 7 < 10®, 1.01 < F < 5, 
0.001 < B < 0.05 G) yield growth timescales of the order of 
a few minutes for the firehose instability. 



4 RESULTS 

We show the results for the maximum growth timescale eval- 
uated using Eq H12|) in Figure 1. We use the bulk Lorentz 
factor F, the random electron Lorentz factor 7 and the value 
of the ambient large scale magnetic field B as our parame- 
ters. We now comment on the range of values we use for the 
parameters employed in the calculations. 



4.1 Range used for 7 

As mentioned earlier, we assume that the TeV radiation is 
produced due to energetic electrons upscattering soft pho- 
tons via the inverse Compton mechanism. Values quoted in 
the literature for the random Lorentz factor 7 used to gener- 
ate TeV emission via the inverse Compton mechanism range 
from « 10* (e.g., Begelman et al 2008) to ^ 10*^ (Ghisellini 
et al. 2009). We use two values for the random Lorentz factor 
of the electrons in our calculations: 7 = 5 x 10* and 10®. 



4.2 Range of magnetic field values 

We consider magnetic field values in the range 0.001 < B < 
0.05 G. This is representative of the values derived for TeV 
blazars. Using a one-zone synchrotron self-Compton model, 
multi-wavelength data fitting requires that the magnetic 
field in the emission regions of TeV blazars MrK 421 and 
MrK 501 be ^ 0.01 G (Abdo et al. 2011a; 2011b). Similar 
considerations yield a value for the magnetic field of 0.03 G 
for the high state of MrK 421 (Shukla et al. 2012). Giroletti 
et al. (2004) obtain B ~ 0.01 G from equipartition consid- 
erations in relation to radio observations of the TeV blazar 
MrK 501. 



4.3 Range used for F 

We next discuss the values we adopt for the bulk Lorentz 
factor F. We use F — 1.01 (corresponding to P||/Pi = 1.32) 
and F = 5 (corresponding to P||/Px = 100). Since the as- 
sumption regarding mild anisotropy (Eq[T|) requires that the 



5 SUMMARY AND DISCUSSION 

We have considered the problem of minute timescale vari- 
ability in TeV blazars. We assume that the TeV radiation 
is produced by highly relativistic electrons (7 ~ 10* — 10®) 
scattering off soft photons via the inverse Compton mecha- 
nism. The jet is beamed towards the observer with a bulk 
Lorentz factor F. These aspects are similar to those con- 
sidered in previous treatments (e.g., Begelman et al. 2008; 
Ghisellini et al. 2009). However, we point out here that bulk 
streaming of the jet plasma results in the parallel pressure 
P|l being > the perpendicular pressure P±. We furthermore 
show that even such a mild pressure anisotropy can give rise 
to a nonresonant, fluid instability called the firehose instabil- 
ity, which can result in a disruption of the large-scale jet over 
timescales of a few minutes. In other words, we have shown 
that the firehose instability is essentially unavoidable in the 
type of situations considered by Begelman et al. (2008) and 
Ghisellini et al. (2009). 

For 7 ranging from 5 x 10* to 10®, F > 1 and 0.001 < 
B < 0.05 G, we find that the firehose instability growth 
timescale is of the order of a few minutes (Figure 1), which 
is in agreement with the observed TeV variability timescales 
for Mrk 501 (Albert et al. 2007) and PKS 2155-304 (Aha- 
ronian et al. 2007). Hence the disruption of the large-scale 
jet features due to the firehose instability provides a natural 
explanation for the observed variability at TeV energies. 
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